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^ | Abstract 

■ It is shown that for moderately hot polar caps (with effective temperature of ~ 10 6 K), the 

efficiency of polar gap acceleration is lower compared to the case in which the polar caps are 
relatively cool and inverse Compton scattering plays no role in controlling the gap. For young 
t-h ■ pulsars with superstrong magnetic fields (> 10 9 T) and hot polar caps (with temperature 

, of > 5 x 10 6 K), because of the energy loss of electrons or positrons due to resonant inverse 

ly-j | Compton scattering in the vicinity of polar caps, pair cascades occur at distances further 

away from the polar cap, and in this case, we have a relatively high acceleration efficiency 
. with ions carrying most of the particle luminosity. 
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,§ ; 1 Introduction 

a 

^ ■ The mechanism for converting rotational energy of pulsars to electromagnetic radiation has long 
• i— i ! been the subject of active research in pulsar physics. The recent observation of seven pulsars by 
the instruments aboard CGRO provides even more challenge to the current pulsar theory. These 

c3 ' pulsars radiate high-energy photons with apparently high efficiency (Ulmer 1994; Thompson et al. 
1994). One possible way to convert the rotational energy to electromagnetic radiation is through 
particle acceleration by the rot at ion- induced electric field in the pulsar's magnetosphere. There 
are several sites for gap acceleration, which include the polar gap (Ruderman & Sutherland 1975; 
Arons & Scharlemann 1979; Michel 1974; Fawley, Arons Sz Scharlemann 1977), the slot gap (Arons 
1983), and the outer gap (Cheng, Ho & Ruderman 1986). In the polar gap model, the acceleration 
region is located near the polar cap while in the outer gap model, the acceleration occurs near the 
light cylinder. In the slot gap model, the acceleration is at the boundary between the open and 
closed field lines. 

There is observational evidence that polar caps of young pulsars can be hot with effective 
temperature T ~ 10 6 K (Ogelman 1991; Greiveldinger et al. 1996). Although the exact value of T 
is rather uncertain, and may depend on the model of the atmosphere above the polar cap (Romani 
1987), the hot polar cap appears to be a plausible consequence of polar cap heating as the result 
of particle acceleration (Cheng & Ruderman 1977; Arons & Scharlemann 1979; Luo 1996). One 
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of the important consequences of hot polar caps is that it provides an alternative mechanism for 
controlling the gap, i.e. mainly through pair cascades initiated by Compton scattered photons 
(Sturner 1995; Luo 1996). Inverse Compton scattering in the polar cap region is strongly modified 
by the magnetic field in that the scattering cross section is enhanced by the cyclotron resonance 
(Herold 1979), i.e. in the electron rest frame, soft photons have energies close to the cyclotron 
energy e B = B/B c in m e c 2 with B c xs A A x 10 9 T the critical field. This process is referred to as 
resonant inverse Compton scattering (RICS). Because of the resonance, the effective cross section 
is greatly enhanced in comparison with the ordinary Compton scattering. For gap acceleration, 
we can define its efficiency as the ratio of the total voltage across the gap to the maximum voltage 
across the polar cap (i.e. for an empty magnetosphere). The efficiency of the gap acceleration is 
controlled by the pair production and thus strongly depends on the effective temperature of the 
polar caps. 

In this paper, we consider the constraint imposed by RICS on the efficiency of polar gap 
acceleration. We assume that free emission of charges from polar caps is allowed, and consider 
both the polar region in which the outflowing charges are electrons (B • fi > 0) and that in 
which the outflowing charges are heavy ions or positrons (B • ft < 0). In the polar gap at which 
B • fl < 0, since the energy loss of ions due to RICS is negligibly small, they can be accelerated by 
the full potential drop across the gap. We derive the condition for the ions to carry most of the 
particle luminosity. Accelerated ions may produce pairs in the anisotropic thermal photon field. 
We calculate numerically the distance from the polar cap at which the ions start to produce pairs 
and compare it to the gap length constrained by pair production by RICS. 



2 Limits on particle energies by resonant inverse Compton 
scattering 

For a hot polar cap which emits soft photons, RICS can constrain the maximum energies of 
electrons (or positrons) accelerated in the polar gap in the following two ways: the energy loss 
competes with acceleration, or by limiting the acceleration length through pair production. For 
a given accelerating potential, <j)(x), the Lorentz factor of an electron (or positron) at a distance, 
xRo, to the polar cap is given by (Dermer 1990; Luo 1996) 

e R a T 6 f 

7(x)=7o + ^ 0(x) -^^io 

where 70 is the initial Lorentz factor, A c = h/m e c is the Compton wavelength, O = 1.686 x 
10~ 4 (T/10 6 K) is the effective temperature (in m e c 2 ), Rq is the star's radius, ot is the Thomson 
cross section, a/ ~ 1/137 is the fine structure constant, £(x') = EB/[y(x')(l — f3(x') cos# max )], 
and ^max is the maximum angle between a photon and the particle motion. For thermal emission 
from the entire neutron star's surface we have # ma x = arccos{[(x + l) 2 — l] 1 / 2 /{x + 1)}, which 
corresponds to propagation of photons at directions tangential to the surface. For a hot polar 
cap, this is valid only for x < 0.5# 2 (where 9 C = (2irRo/ Pc) 1 ^ 2 and P is the pulsar period); for 
x > 0.59 2 , the angle is given by 9 max ~ arccos[a;/(6' 2 + x 2 ) 1 ^ 2 ], corresponding to photons from the 
polar cap rim on which the field lines have colatitude angle 9 C relative to the magnetic pole. The 
energy loss due to RICS is described by the third term on the right hand side of equation (1). 

The maximum Lorentz factor corresponds to the value at x — x max which is the gap length 
determined by the pair production, e.g. due to RICS or due to thermal photons interacting with 
heavy ions, whichever comes first. Equation (1) can be solved numerically with the boundary 
condition 7(0) = 70 . 
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The model for the accelerating potential strongly depends on the physical conditions on the 
polar cap. If polar caps are sufficiently hot, we may have free emission of charges. The induced 
electric field causes charge flow, which tends to screen out the field. However, the charge screening 
is always insufficient, e.g. due to field line curvature (Arons & Scharlemann 1979), and hence 
charge depletion occurs (the charge density deviates from Goldreich- Julian (G-J) density). An 
electric potential exists in the charge-depletion region (often called the gap). As an example, we 
consider the potential (e.g. Arons & Scharlemann 1979; Arons 1983) 



<f)(x) « 0.3759*cBRo ( — 



2i 



g(x) sin^ sini, (2) 



where 9 is the angle of the field lines relative to the magnetic pole, ip is the azimuthal angle of 
the open field line flux tube, x is the distance to the polar cap (in Rq), i is the angle between 
the rotation axis and the magnetic pole, and g(x) is given by g(x) = x 2 for x < 0.59 c , g(x) = 
6 c \(x + l) l l 2 — 1] for x > 0.5# c . According to Arons (1983), the angle ip is in the ranges, — ir < ip < 
for % > ir/2 (if the outflowing charges are ions), and < ip < ir for % < tt/2 (if the outflowing charges 
are electrons). The potential given by (2) is due to charge density decrease along field lines (i.e. it 
becomes less than the G-J density) because of field line curvature. Plots of the Lorentz factor as a 
function of distance from the polar cap are shown in Figure 1. We assume 7 = 10 at xR = 1 cm. 
As shown in the figure, energy loss due to RICS can compete with acceleration only for high polar 
cap temperature T > 5 x 10 6 K. For T ps 10 6 K, the particle energy is constrained mainly through 
pair production by RICS. 

There are other models of the space-charge-limited potential such as the model discussed by 
Michel (1974), also by Fawley, Arons & Scharlemann (1977) and Cheng & Ruderman (1977), which 
do not rely on field line curvature and in general gives a lower voltage. The potential discussed 
by Michel (1974) arises from the deviation of the charge density from the G-J value because of 
particle inertia, i.e. the velocity v(x) and the magnetic field B(x) have different scaling with x 
such that the charge density cannot be maintained at the G-J density all the way along open field 
lines. The corresponding potential can be described by (Fawley, Arons & Scharlemann 1977) 

<t>(x)=2xcBR [—) , (3) 

with Q = 2n/P, Qb — ZeB/m and x < 9 C . For x > 9 C , the electric field drops off faster than 
(xRq)~ 2 . The plot of the Lorentz factor as a function of xRq is shown in Figure 2 for T = 1.5 x 10 6 K 
and T = 5 x 10 6 K. The electric potential for x > 9 C is assumed to be constant, 4>(x) = <f){0 c ). For 
T = 5 x 10 6 K, the energy loss is so severe that it continues beyond the distance x > 9 C , resulting 
in a sudden decrease in 7 at x = 9 C . As in Figure 1, for moderately hot polar caps (T 10 6 K), 
the particle energy is constrained by electron-positron cascades started by RICS. 

The main mechanism for a high energy photon (emitted by an electron or positron either 
through curvature or inverse Compton scattering) to produce a pair is single photon decay in 
a strong magnetic field. To produce a pair, the photon must satisfy the threshold condition 
0.5esin6> 7 > 1 where e is the photon energy (in m e c 2 ), 9 1 is the angle between the photon prop- 
agation and the magnetic field. When this condition is satisfied, the number of pairs produced 
depends on the opacity, r, which is an integration of the absorption coefficient (for a photon being 
converted into a pair) along the distance that photon travels. We assume that a photon produces 
at least one pair if r > 1. With this condition, we can write down the threshold condition for 
producing at least one pair through RICS (Luo 1996) as 

, n 7 15\/ B \- 2 f P \V2(i + x )6 

^ >134 (toA)(37To^) (an) (4) 
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Figure 1: Lorentz factor versus distance to the polar cap. The electrons (or positrons) are 
accelerated by the potential (2) with i = tt/A (or 3tt/4), B = 3 x 10 9 T, P = 0.15s. The 
dashed curve is obtained without RICS. The solid curves (from top to bottom) correspond to 
T = 1.5 x 10 6 K, 5 x 10 6 K, and 10 7 K. We assume that thermal emission is from polar cap. The 
dotted curve is the threshold condition for pair production through RICS (i.e. pair production 
can occur only above the curve). 
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Figure 2: As Figure 1 but the particles are accelerated by the potential (3) with % = (i = tt for 
outflowing positrons or ions), P = 0.15 s, B = 3 x 10 9 T, T = 1.5 x 10 6 K (the upper solid curve), 
5 x 10 6 K (the lower solid curve). The result is similar to that derived by Sturner (1995). The two 
dashed plots (upper and lower) correspond respectively to acceleration of electron (positron) and 
ion without RICS. The dotted curve is the threshold condition for pair production through RICS. 
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Figure 3: Pair production due to ion colliding with soft photon from polar cap emission above the 
polar cap with ft ■ B < 0. Beside ions, the outflowing charges may include positrons. 



The right-hand side of (4) strongly depends on B. If the above condition is satisfied, pair cascade 
can be initiated by RICS. The threshold condition is represented by the dotted curves in Figures 
1 and 2. Note that for each curve of log 7 vs. log(xi?o/cm), the location of an intersect (with the 
dotted curves) point defines a gap length, represented by Xq 10S . (If the curve has two intersecting 
points, the one that is the nearest to the polar cap is relevant.) 

3 Pair production by photons in the field of the nucleus 

Electron positron pair cascades can be started by accelerated ions in thermal photon fields. A 
photon can decay into an electron positron pair in the Coloumb field of a nucleus when the center 
of momentum frame energy exceeds the rest mass of the nucleus plus two electrons, that is 

s > (m N c 2 + 2m e c 2 ) 2 , (5) 

where m N is the mass of the nucleus. This condition can be satisfied for soft photons from the 
polar cap in the presence of accelerated ions. We can estimate the mean free path of an ion before 
producing an average one pair. If we believe one pair per primary particle is enough to short 
out the electric field, this mean free path length, represented by x™ n , can be identified as the gap 
length (Figures 3 and 4). 

Numerical evaluation of the mean free path for an accelerated ion traversing an isotropic 
thermal photon field was given by Protheroe (1984), and Protheroe & Johnson (1996). The 
injection rate of pairs for a given energy of iron nuclei was numerically calculated by Bednarek & 
Karakula (1995) for an anisotropic photon field from the polar cap. We calculate number of pairs 
produced by iron nuclei accelerated along open field lines by the potentials given by Eqs. (2) and 
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Figure 4: Plots of cumulative number of pairs produced per ion as a function of distance. Solid, 
dashed and dotted curves correspond respectively to T = 5 x 10 6 K for the whole star's surface, 
T = 10 7 K for the polar cap (but the star's surface has T = 10 6 K), and T = 5 x 10 6 K for the 
polar cap (but the star's surface has T = 5 x 10 5 K). The thick and thin curves are calculated 
using the potentials (2) and (3), respectively. We assume P = 0.15 s, B = 3 x 10 9 T. 
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(3), and the mean free path corresponds to the average distance at which one pair is produced per 
ion nucleus. The mean cumulative number of pairs produced in the radiation field of the polar cap 
(or whole neutron star) at temperature T by a nucleus of energy E traveling along the neutron 
star magnetic axis from the pole to distance xRo above the polar cap surface is given by 

T? ex roc /-cos[0 m i n (£,e)] 

N e+e -(E)^^- / n(e) (1 - /3cos 0)a(s) dcos0 de dx, (6) 

2 Jo Jo Jcos[O max {x)] 

where n(e) is the differential photon number density of photons of energy e in black body radiation 
at temperature T, and <j(s) is the total cross section for pair production (Protheroe 1997) by a 
the nucleus at center of momentum (CM) frame energy squared, s, given by 

s = m 2 N c 4 + 2eE{\-(3 cos 6) (7) 

where 9 is the angle between the directions of the nucleus and photon. Note that 6 min (E,e) is 
determined by the threshold condition (5). Plots of accumulative number of pairs produced per 
ion through interaction with thermal photons as a function of distance from the polar cap are 
shown in Figure 4. for P = 0.15 s and B = 3 x 10 9 T. Compared to RICS, in this case we have 
Xq ics < x on . Although we use a high B, the process strongly depends on P and T. 



4 Acceleration efficiency 

The efficiency of the polar gap is described by the ratio of the potential drop across the gap, 4>(x ), 
to the maximum potential across the polar cap with an empty magnetosphere, max , viz. 

v = ^) , (8) 

where max = 0.5# c 4 aB-Ro ~6.6 x 10 12 V (1 s/P) 2 (£?/10 8 T), 6 C = (flR /c) 1/2 = 1.45 x HT 2 (ls/P) 1 / 2 , 
and x is the length of the polar gap in units of R . Let L p and L R be the accelerated particle 
and spin-down luminosities, respectively. If the energy loss of accelerated particles in the gap is 
not important, we have r\ L p /L r since we can write L p ps Ze(ft(xo)N GJ and L R ^ Ze0 max iV GJ 
where N GJ is the injection rate calculated from the G-J density n GJ ~ 2eoBQ/(Ze). Note that 
the efficiency defined by (8) is similar to that used by Arons (1996). For the presently available 
models for accelerating potentials, we always have rj < 1. 

Using (8), we may estimate the efficiency for a given model potential. For the potential 
described by Eq. (2) and assuming x ^> 0.5# c , we have the maximum efficiency 

^max ~ 0-296 c x^ 2 sinip sfiu. (9) 

Thus, even when we use Xq ~ Rl/Rq = I/O 2 , the maximum efficiency is ^ max ~ 0.29. When a pair 
cascade occurs, we usually have xo <^ 1/9 2 (for the polar gap) and hence i] i] max . 
One may estimate the efficiency for the potential given by Eq. (3) from 

(A\ l / 2 (x\ /10 8 T\V2/px 
77 = 0.069 - - — 77 , (10) 



z) \ej\ b ) Vis 

with x < min{# c , x$}. Although we have a higher rj for the Arons & Scharlemann model, the 
acceleration by the potential given by Eq. (3) is the more effective for i m 0. 

The evaluation of xq depends on the specific mechanism for initiating a pair cascade and 
whether the pair plasma produced through the cascade is dense enough to short out part of the 
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electric field. We consider three mechanisms for starting a pair cascade: RICS, curvature radiation, 
and inelastic scattering of ions by thermal photons from the polar cap. The corresponding lengths 
are Xq ics , x™ rv and x 1 ™. In general, we have Xq ics < XQ mv (e.g. Luo 1996). For moderately 
hot polar caps with effective temperature ~ 10 6 K, the energy loss due to RICS is not important 
compared to acceleration but the photons produced through RICS can start a pair cascade at the 
distance less than Xg ICS < Xq utv , x 1 ™. Thus, the effect of RICS is to reduce the gap acceleration 
efficiency For space-charge-limited flow, the composition of outflowing charges at the poles with 
f2 • B > and Q ■ B < can be different. For Q • B > 0, the primary particles consist mainly of 
electrons, and for O • B < 0, the main components are heavy ions or positrons. In the ion zone, 
the gap length is constrained by the pair production by positrons through RICS. The possible 
source of positrons was discussed by Cheng & Ruderman (1977). For moderately hot polar caps, 
the gaps at both types of pole have a similar efficiency. 

For sufficiently hot polar caps and a superstrong magnetic field (~ 10 9 T), the energy loss 
of electrons or positrons due to RICS can be important, and may prevent them starting a pair 
cascade in the region close to the polar cap. Since the cross section is oc (Z 4 / ' A 2 ){m e /m p ) 2 ps 1CT 4 
for iron nuclei (Z/A 1/2), the energy loss due to RICS for ions is negligibly small, and they can 
be continuously accelerated until distance x < min{xo ICS , Xq 011 }. Thus, the acceleration efficiency 
increases significantly. As an example, for the potential (2) with B = 3 x 10 13 and T = 1.5 x 10 6 K, 
the gap length is about Xq ics Rq ~ 3 x 10 3 cm and the maximum energy of electrons (or positrons) 
is 10 4 m e c 2 . This gives the particle luminosity L p ps 10 4 m e c 2 iV GJ where N GJ is the injection rate 
of primary electrons (or positrons). For hot polar caps with T = 10 7 K, from Figures 1, we 
obtain Xq ics Ro ~ 1.6 x 10 4 cm. From Figure 4, pair production by ions interacting with thermal 
photons is important only for the case where the whole star's surface has T = 5 x 10 6 K (the 
thick solid curve), and this occurs at the distance a;o° n i?o ~ 10 5 cm. Thus, the gap length is 
controlled by RICS. From Figure 1, we have 7(xo ICS ) « 1.6 x 10 5 for positrons (or electrons) and 
7(xo ICS ) ~ (Zm e / Am p ) x 10 6 for ions; positrons (or electrons) and ions can be accelerated to 
the maximum energy 1.6 x 10 5 m e c 2 and 10 6 Zm e c 2 , respectively. Hence, we have much a higher 
particle luminosity L p ps 10 6 m e c 2 A^ G J and most of L p are carried by ions. 

5 Application to young pulsars 

From observations, several pulsars may have polar caps with effective temperature > 10 6 K though 
the exact value of the temperatures is rather uncertain (for a review, see e.g., Ogelman 1991). 
For moderately hot polar caps and typical B fields (B ~ 10 s T), the energy loss rate is not rapid 
enough to compete with the acceleration, and RICS controls the polar gap only through pair 
production. This in general reduces the gap length, and hence reduces the acceleration efficiency 
by at least one order magnitude, imposing a severe constraint on the energetics of the polar gap. 

For young pulsars with superstrong magnetic fields > 10 9 T, the gap acceleration efficiency 
can increase provided that the polar cap (O • B < 0) is hot (T > 5 x 10 6 K) and for electrons 
(or positrons), the energy loss due to RICS can dominate over the acceleration. One of CGRO 
pulsars, PSR 1509-58, has magnetic field > 1.5 x 10 9 T (the magnetic field at the pole can be 
twice this value). Since neutron stars are hot at birth (after supernova explosion), their polar 
caps can be hot as result of self-sustained polar cap heating at very young age (e.g. Luo 1996), 
e.g within hundred years or so after birth. These high field pulsars, assuming they were born 
with high magnetic fields, must have undergone a phase when thermal emission had an important 
effect on pulsar electrodynamics. 

We suggest that for high field pulsars there can be four phases during the their lifetime as 
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discussed below, (i) At very young age (within a few years after supernova explosion) when 
the surface or polar cap temperature is high, energy loss due to RICS is so severe that a pair 
cascade due to electrons or positrons cannot occur in the region close to the polar cap. In the 
ion zone (ft • B < 0), accelerated ions may carry most of the particle luminosity. For heavy 
ions, photodisintegration may occur in the thermal photon field from the polar cap. At this 
stage, since pulsars also have a very large spin-down luminosity, they are a potential source of 
high energy hadrons and electromagnetic radiation, (ii) As pulsars cool down, the energy loss 
rate due to RICS becomes less than acceleration rate. However, RICS may still be the dominant 
mechanism for pair production to constrain the gap length (Luo 1996). The high field pulsar, 
PSR 1509-58, may currently be at this stage since the polar cap temperature inferred from the 
observation is about 1.5 x 10 6 K (e.g. Kawai et al. 1991). Its gap length is shorter than that 
controlled by pair production by either curvature radiation or relativistic ions in the photon field. 
Thus, the acceleration efficiency is less than in phase (i). (iii) As the surface temperature or 
polar cap temperature becomes much lower, pair production by RICS is not important. The 
main mechanism for initiating pair cascades is curvature radiation. As the period increases, the 
gap length increases and hence the efficiency also increases. However, at this stage, the pulsar 
spin-down luminosity is lower than at their early young age. (iv) The period reaches the critical 
period at which pair production due to curvature radiation is no longer operative. Although the 
gap acceleration efficiency is now high (no pairs to screen the electric field), the luminosity is very 
low. 

In summary, for very young pulsars with superstrong magnetic fields (B > 10 9 T) and hot polar 
caps (> 5 x 10 6 K), the severe energy loss due to resonant inverse Compton scattering can prevent 
electrons or positrons starting a pair cascade in the region close to the polar cap but ions can 
be accelerated by the full potential drop across the gap (provided that free emission of ions can 
occur). We then have the ions carrying most of the particle luminosity. Thus, very young pulsars 
with superstrong B field can be a potential source of high-energy hadrons and gamma-rays. In 
contrast, for pulsars with moderately hot polar caps (T ~ 10 6 K), the effect of resonant inverse 
Compton scattering is to reduce the gap acceleration efficiency Pair production by accelerated 
ions interacting with thermal photons strongly depends on P and T, and in general we have 
Xq 10S < and hence the polar gap length is controlled by RICS. 
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